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_ Abstract 

D ' Unexpected features in the energy spectra of cosmic rays electrons and positrons have been 

recently observed by PAMELA and Fcrmi-LAT satellite experiments, opening to the exciting pos- 

O^l , sibility of an indirect manifestation of new physics. A TeV-scale fourth lepton family is a natural 

' extension of the Standard Model Icptonic sector (also linked to the hierarchy problem in Minimal 

^Kf^ ' Walking Technicolor models). The heavy Majorana neutrino of this setup mixes with Standard 

' Model charged leptons through a weak charged current interaction. Here, we first study analyti- 

, cally the energy spectrum of the electrons and positrons originated in the heavy Majorana neutrino 

decay modes, also including polarization effects. We then compare the prediction of this model 

with the experimental data, exploiting both the standard direct method and our recently proposed 
- 1 — I 

I Sum Rules method. We find that the decay modes involving the tau and/or the muon charged 

, leptons as primary decay products fit well the PAMELA and Fermi-LAT lepton excesses while 

there is tension with respect to the antiproton to proton fraction constrained by PAMELA. 
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I. INTRODUCTION 



A hardening in the electron and positron fluxes has been observed in Cosmic Rays 
(CRs) at energies between 10 GeV and 1 TeV. More precisely, PAMELA [l| measured 
a rise in the positron fraction above 10 GeV. Fermi-LAT 0] measured the sum of the 
positron and electron fluxes, detecting a feature with spectral index of about -3 above 100 
GeV, followed by a decrease at 500 GeV. The latter results strengthened previous hints 
by ATIC 0] and PPB-BETS 0] balloon experiments and by the ground-based telescope 
HESS fl. 

When compared with the predictions of current astrophysical models, the results of 
PAMELA and Fermi-LAT independently suggest the presence of an excess of electrons 
and/or positrons. The origin of these excesses is unknown and many explanations have 
been put forward: from misunderstood astrophysical effects to new physics indirect effects. 
Notice that the lack of any excess in antiproton CRs found by PAMELA is a crucial 
element in understanding the origin the electron and positron excesses. Interpreting these 
excesses as signatures of a dark matter component, generically an annihilating or decaying 
WIMP, is for sure exciting; see e.g. [3| for a review. 

While the interpretation in terms of dark matter annihilations often leads to an un- 
observed excess of gamma and radio photons, the interpretation in terms of dark matter 
decays is compatible with photon observations 0, 0] , even though some channels now start 
to show some tension [10l4l3l|. A discrimination strategy was proposed in 

In this work, we want to investigate a simple extension of the Standard Model (SM): 
a heavy Majorana neutrino which couples to SM particles only via charged current 
interaction: 

cf^ = Q w-e^PLN + h.c. (1) 

where Cg, with i = e, /i,r, is a numerical factor parametrizing the strength of the mixing 
between the heavy neutrino and the SM charged lepton i and Pl is the L-chirality 
projector. The interaction above naturally arises extending the SM leptonic sector by 
introducing a fourth lepton family, as happens e.g. in Minimal Walking Technicolor models 



mm]- 

We are interested in a nearly stable and TeV-scale heavy neutrino, which could con- 
stitute a fraction of the dark matter density, xjy = pn/pdm- Via its semileptonic decays 
N —7- i'^W^ - whose branching ratio (BR) into the specific flavor £ is proportional to C| 
- the heavy Majorana neutrino could also be a source of electron and positron CRs. Our 
goal is to study whether the electron and positron excesses observed by PAMELA and 
Fermi-LAT can be attributed to such decaying heavy Majorana neutrino. Notice that in 
our model the fluxes of electrons and positrons coming from the decay chains are equal. 

We focus on a very specific model: we are thus able to discard the model if it turns out 
not to fit the data, or to strongly constrain its parameter space. Instead of using numeri- 
cal codes 11, 22 -25 1, for each lepton flavor £ here we calculate analytically the energy 
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spectra of the electrons and positrons produced in — t- iW and the subsequent £ and W 
decay chains, including polarization effects. This allows to have a deeper understanding of 
the physical results and, as we are going to discuss, to find features that are not properly 
accounted for in some numerical codes. We then consider the propagation of electrons and 
positrons according to the most popular models and compare our predictions with the ex- 
perimental results of PAMELA and Fermi-LAT, exploiting both the more straightforward 
Sum Rules (SR) method proposed in ref.[26| and the longer direct comparison method. 

Considering in turn each channel N — )• iW (with i = e, fi, r) and its subsequent decays, 
we find that: the decays into tW fit the experimental data extremely well for M^r ~ 3 
TeV and , a good fit is obtained for the decays into iJ,W for the same model 

parameters; while it is impossible to reproduce the data if decays mostly into eW. 

As for the CR antiproton contribution from the heavy Majorana neutrino decays, in 
ref . 1241] it was estimated to be compatible with the PAMELA experimental data. In ref. 



there was some tension with the data, but a fit to the ATIC data were used there. 
We have redone the analysis for the antiproton flux as well as the antiproton to proton 
fraction for the processes relevant here and compared to CAPRICE 27] and PAMELA 0] 
data. We find that it is possible to accommodate the CAPRICE data while we observe a 
tension with the PAMELA data. 

We conclude that a 3- TeV fourth family Majorana neutrino with dominant BR into 
tW and/or fiW with coupling x^^'^C^/^ ~ 10^^'^ represents a plausible and even quite 
conservative explanation of the electron and positron CRs excesses. 

Our findings confirm previous numerical studies for the e and /i channel, but disagree 
for the T channel, which was considered to be disfavored [2j] or to provide a worse fit to 



the data as compared to the fi channel |lll ]. These analysis were based on the Monte Carlo 
simulation program PYTHIA which (as well as HERWIG) treats leptons and vector bosons 
as unpolarized. To understand the origin of this discrepancy, we carried out our analysis 
also neglecting polarization effects. We find that, although going in the right direction, 
polarization effects cannot alone explain the discrepancy between our analytically-based 
results and the numerically-based results of previous analyses. A comparison with ref. 0] 
is not strictly possible, since polarization effects were neglected and especially the ATIC 
data were used for the fits, not the subsequent ones by Fermi-LAT Q]. On the other 
hand, in the more recent analysis of ref.|25l] including polarization and even electroweak 



emission 



281 ] ■ the semileptonic decay channels were not considered. 



The paper is organized as follows. In section II the model setup is introduced. In 
section III the analytical calculation of the energy spectrum of electrons and positrons 
from A'^ decays is carried out, studying separately each of the three semileptonic channels. 
Section IV is devoted to the comparison of the model prediction with the experimental 
data, bot with the SR and direct comparison methods. Section V presents the estimate 
for the antiproton flux and antiproton-proton ratio. We draw our conclusions in section 
VI. 
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II. THE FOURTH LEPTON FAMILY SETUP 



In the SM, the three lepton families, £ = e, /x, r, belong to the following representations 
of the gauge group SU{3)c x SU{2)l x U{1)y. 



Li = {yiL (1,2,-1/2) 



(1,1,-1) 



(2) 



where the chirality projectors Pl = (1— 75)72 and Pr = (1+75)72 have been introduced. It 
has been observed experimentally that at least two of the SM neutrinos have a small mass, 
not larger than the eV-scale [29|. In the following, we account for light neutrino masses 
and mixings via an effective Majorana mass term, i.e. by adding to the SM Lagrangian a 
dimension-5 non-renormalizable operator (which could arise for instance through a seesaw 
mechanism) . 

Our aim here is to investigate the CR positron and electron fluxes originated via the 
decay of a possible TeV-scale fourth lepton family, for which we introduce the ("-flavor. 
This additional family is composed by a lepton doublet, a charged lepton singlet and a 
gauge singlet: 



Lr 



(^C 



(1,2,-172) 



:i,i,-i) 



(3) 



We keep our scenario as general as possible, but we stress that it would arise naturally 
in Minimal Walking Technicolor models [la], with interesting collider signatures 19| - see 



also 



20( 1 and, for a complete review, ref. [211]. 

The (-charged lepton, C = Ci + C-Ri h^-s a Dirac mass term like the other three charged 
leptons of the SM, but large enough to avoid conflict with the experimental limits. We 
work in the basis in which the 4x4 charged lepton mass matrix is diagonal. The mass 
eigenstates of the two electrically neutral Weyl fermions instead correspond to two Majo- 
rana neutrinos, A^i,25 for which we are free to choose the ordering Mjv^ < M^Vj. In such a 
framework, assuming the C, lepton to be heavy enough with respect to A'^i, it turns out that 
the latter is stable and can constitute - at least part of - the dark matter in our galaxy. 

We consider the possibility that the new heavy leptons mix with the SM leptons. 
For clarity of presentation, we assume that the heavy neutrinos mix only with one SM 
neutrino of flavor i {i = e,n,T). This is not a restrictive assumption: as demonstrated in 
appendix lAl the expressions that we are going to derive apply also to the case in which 
the heavy neutrinos mix simultaneously with all three light lepton families. The entries of 
the mass matrix are: 



^ 0{eV) 0{eV) \ / 



0{eV) 0{eV) mo 
mi ruD ruR j 



h.c. 



(4) 



The measured values of the light neutrino masses suggest that the entries of the upper 
2x2 block have to be of 0{eV). The mass scale of niR could be as large as the cutoff of the 
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theory, since it is not protected by any symmetry (we allow violation of the lepton flavour 
number by 2 units). The elements m£ and uid are expected to be at most of the order 
of the electroweak scale. Given such a hierarchical structure and up to small corrections 
of 0{eV/MNi,N2) ^ 10^^^, one obtains the following form for the unitary matrix which 
diagonalises eq.([l|): 



\ J 



V 



Nil 

\N2L j 



V 



( 



\ 



cos fli I COS c/ sm f 
- sin Qn i cos 6 cos ( 
—isinO 



sm ti sm f 
sin 9 cos / 
cos 9 



\ 



(5) 



-^£,1,2 are the new Majorana mass eigenstates (actually A'^ is a flavor eigenstate when 
including all three SM leptons, as discussed in the appendix) and 



tan 9e 



me 
rriD 



tan 261 



rriR 



(6) 



The light neutrino Ni has a mass of 0{eV). Up to corrections of 0{eV), the heavy 
neutrinos Ni^2 have positive masses given by: 



1 + 4 



1 



m 



'1 + 4- 



mjj + 



+ 1 



R 



m 



(7) 



R 



The smaller is rriji (namely the more 9 approaches vr/4), the more the neutrinos A^i and 
become the two Weyl components of a Dirac state. Notice that the mixing between 

the heavy neutrino states A^i^2 is controlled by the mixing angle 9£, which is constrained 

experimentally to be small (more on this later). 

The neutral current interactions of the neutrinos are 



'-N 



a 



4 cos 9. 



-Z^{Na^-f5Ni + cos^ 9 A^i7^75A^i + sin^ 9 iV27^75A^2 - i sm{29) A^27^A^i) , 

(8) 

while their charged current interactions with the light and heavy charged leptons are 
respectively 



rCC 



g_ 

^/2 



W"Cl7^(- sin 9iNiL + icos9 cos 9iNiL + sin 9 cos 9iN2L) + h.c. (9) 



and 



Cf^ = — ^W-£Ll^{cos9iNi,L + icos9sm9iNiL + sm9sm9eN2L) + h.c. . (10) 
v2 

Notice that the neutrino neutral current remains flavor diagonal at tree-level (the neutral 
current is not flavor diagonal in models with TeV scale right handed neutrinos involved in 
the see-saw mechanism for the light SM neutrino masses, see e.g. [30(]), hence the heavy 
neutrinos couple to the SM ones only through the charged current interactions at this 
order. This is a distinctive feature of our fourth lepton family. 
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Indeed, also the Yukawa interactions do not mix the heavy Majorana neutrinos with 
the SM one: 

= — — — - — (iViA'i + N2N2) - ? cos(26')A'i75A'2 H . (11) 



2Mw 

In this work we are then concerned with a heavy Majorana neutrino A'^i, which only 
couples to the SM charged leptons, and with a small mixing angle. If such mixing is small 
enough, A^i can be considered as nearly stable on cosmological time scales. We assume 
> Ml + M\Y, so that the charged lepton C already disappeared as a dark matter 
candidate because of its decays into WNi. 

As for the possible hierarchies between the heavy neutrinos, two main situations arise: 
if M7V2 > Mtvj + Mz, also already disappeared; if Mjy^ < M^r^ + Mz, N2 is nearly 
stable like A'^i (in the limit of degenerate masses they form a Dirac neutrino). From now 
on, we consider the first scenario. 



A. Laboratory constraints 



From collider experiments there are both direct and indirect constraints on heavy 
leptons. We refer to |l9l | for an updated analysis (together with some useful references) 
and here we just summarize the results relevant for the present analysis. 

As stressed, we are concerned with a heavy neutrino A^i which is nearly stable on 
cosmological time scales. For this setup, a lower bound on the heavy neutrino mass Mtvj 
can be extracted from the constraints on the effective number of neutrinos involved in Z 
decay: M^r^ > 44(45.4) GeV at 2a for a Majorana (Dirac) nearly stable heavy neutrino. 
Note however that in the present analysis we are interested in much bigger values of M^-^ , 
as will be discussed in the following sections. 

As for the heavy charged lepton the lower limit on its mass is 102.6 GeV if it is stable. 
This limit is only slightly weakened in the case of a decaying heavy charged lepton, as in 
our setup. 

The mixing angle Oi of the heavy leptons with the SM charge lepton i is strongly 
constrained by lepton universality tests: sin^ 6*^ < 0.012(0.15) at 3a for £ = e, fi {£ = r). 
However, as we are going to discuss in the following, much smaller values of 9e are required 
in order to have Ni as a nearly stable heavy neutrino. Hence, at first order in 9^ the mixing 
matrix V of eq.dS]) becomes: 

/ 1 icOS^sin^^" «inflRinfl„\ 



V 



sm ( 



I cos V COS I 



sm sm I 



sm a cos ( 



(12) 



y —isinO cos 9 j 

We emphasize that this structure arises also in the case of three simultaneous mixings 
with the three SM lepton families, as discussed in the appendix. 
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B. Comological constraints 



It is well known that a stable or nearly stable fourth generation Dirac neutrino is 
excluded as the main component of dark matter. Its vector coupling to the Z results in an 
unsuppressed annihilation cross section into fermion pairs at low momenta, as well as in a 
large elastic scattering cross section with nucleons. The latter strongly constrains the relic 
density of Dirac neutrinos to be several orders of magnitude smaller than the observed 
dark matter density. 

Some of these problems can be avoided when the Dirac pair is split using a Majorana 
mass term as in our setup. The axial coupling of the Z boson with Ni {i = 1, 2) suppresses 
its annihilation cross section into fermion pairs at low momenta, as well as the elastic 
scattering cross section with nucleons. For the vector coupling between Ni,N2 and Z, 
such suppression is not at work. By requiring that M2 — Mi > 0(100) GeV (the typical 
momentum transferred in dark matter-nucleon scatterings) , large inelastic scattering rates 
are avoided 



3ii 



However, even in this case the thermal relic density of A^i, pNi , turns out be significantly 



smaller the dark matter one, pdm- see e.g. [33] for an updated analysis, where it is shown 



that the A^i thermal relic density can reach up to 20% of the observed dark matter density, 
but is below 10% for most of the parameter space. 

We then account for the possibility that pN-^ is smaller than the dark matter one, poM- 
It is however reasonable to assume that they have the same (spherically symmetric) profile: 

PNi{r) = PDAlir) . (13) 

In the following we will consider free parameter and adopt for definiteness the 



Navarro Prenk White (NFW) 33|] profile for the dark matter halo of the Milky Way: 



pDAiir) = po — (14) 

(r/r,)7[l + (r/re)"] — 

where a = 1, /3 = 3, 7 = 1, rc = 20 kpc and po = .26GeV/cm^. The latter normalization 
is chosen in order to have pDM{fsun) = 0.30GeV/cm^, where rgun = 8-5 kpc. 



C. Contribution to electron and positron CRs from A^i decay 



We now turn to the study of the production and propagation of the electrons and 
positrons in CRs, considering the contribution from the decay of the heavy Majorana 
neutrino A^i. 

The number densities of electrons and positrons per unit energy satisfy the transport 
equation 0, [34 - 38 ] : 



J^^{E,±,r) = K{E,±)V^ f,±{E,±,r) + —^[h{E,±)f,±{E,±,r)] + Q,±{E,±,r) , (15) 
ot oE^± 
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where E^± denotes the electron or positron energy. The first term describes the propaga- 
tion through the galactic magnetic field and the diffusion constant is given by K{E^±) = 
Ko{Eg± /GeY)^ . For definiteness, in what follows we adopt the MED-model parameters for 
the propagation, namely Ko = 0.0112kpc^/Myr and 6 = 0.70 Positrons lose energy 
through synchrotron radiation and inverse Compton scattering on the cosmic microwave 
background radiation and on the galactic starlight at a rate h{E^±) = E'^±/{Gq\te) where 
TE = 10^^ sec. Finally, the last term in the equation above represents the source term. 

The assumption that positrons and electrons are today in equilibrium implies that we 
can impose the steady state condition df^±{E^±,r)/dt = 0. The galaxy is described as a 
cylinder of radius R and half-thickness L, for which we assume the values R = 20 kpc and 
-L = 4 kpc respectively. We impose that the number density of positrons and electrons 
vanishes at the surface of the cylinder. 

The source term must include all possible mechanisms through which are produced, 
from astrophysics to possible new physics. As for astrophysics, most CR electrons are 
likely to come from supernovae remnants, while positrons are mainly produced in hadronic 
processes when CR protons collide with intergalactic hydrogen. Clearly, the linearity of 
the diffusion equation guarantees that its general solution can be expressed as a sum of 
the solutions obtained by considering separately each possible source mechanisms. 

As already discussed, we assume that M2,M^ >> Mi, so that and C have already 
decayed into A^i. On the contrary, A^i could be nearly stable since it mixes only with light 
leptons (see eq. ffTO]) : 

where the coefficient Ce = cos9s'm9i has to be small because so is the mixing angle 6i. 

Clearly, for electrons and positrons the source term of eq. (fT5]) associated to the decay 
of the heavy Majorana neutrino A'^i is: 

Qe± {E^± , r = — — (17 

Mn^tn^ dE^± 

where dN/dE^± represents the number of per unit energy that are produced in the 
decay of A^i, t^^ is the lifetime of Ni and pjy-^ is its present relic density. We can split the 
number of per unit energy into 

where dN^^^\E^±) / dE^± refers to the specific channel A'^i — )■ iW (whose vertex is pro- 
portional to ) and the subsequent decays of I and VF, see figlSl Introducing the partial 
lifetime = tm^/ BR{Ni — t- IW)^ the source term can be simply written as a sum of the 
contributions of the three light lepton flavors: 

Q..{E,.,r) = Y.Q'AE^-^r-) , QU{E,.,r) = ^f"] dNi^J\E,.^ . (19) 



= -iJi=c,W-hl''NiL + h.c. , (16) 
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Since at tree level the Feynman amplitudes for the two charge conjugated decays of Ni are 
equal, the energy spectra of positrons and electrons are also equal: dN^^^\E^+)/dE^+ = 
dN^^^\E^-)/dE^-. From now on we will omit the apex it, unless where necessary. 

Clearly, we are interested in fe{Ee,r, z) for z = and r = rgun = 8.5 kpc. The solution 
of the transport equation at the Solar System, can be formally expressed by 

fe{Ee) = feiEe) , /f (i?e) = dE'^ G{E^, K) ^N^'J^iK) (2O) 

where xni = PNi/pdm and G{Ee,E'^) is a Green function which contains all the astro- 
physical dependencies and whose explicit form can be found in 37[ • Notice that the Green 
function is not very sensitive to the choice of the dark matter halo profile, since the Earth 
receives only electrons and positrons created within a few kpc from the Sun, where the 
different halo profiles are very similar. 



In the present analysis we adopt the approximation of 221 ]: 



G{E, E') = l^e'^+^(^'-'-^"")0(i?' -E)—^ (21) 
E^ cm.-^ 

where the coefficients a and h are the appropriate ones for the MED-model and NFW 
profile: a = —1.0203, b = —1.4493. It was found that this approximation works better 
than (15 — 20)% over the whole range of energies 
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III. ELECTRONS AND POSITRONS FROM HEAVY MAJORANA NEUTRINO 
DECAY 

We now turn to calculate the A'^i lifetime and the number of electrons and positrons 
produced in its decay. 



A. Lifetime of A^i 

The decay width for A''i is 

r^vi = 5^ , = r[iVi ^ w+r] + t[Ni ^ w~t] , (22) 

where the two charge conjugated width are equal (at tree level) and, in the approximation 
<C Mw, are given by 

r[iV. W-i-] = f 1 - ^ Vl + ^ - . (23) 
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500 1000 1500 2000 2500 3000 



Mni [GeV] 

FIG. 1: Partial lifetime of A^i as a function of its mass MjVj and for selected values of the combi- 
nation of mixing angles — cos sin 0^. 

The branching ratio shave a simple expression in terms of the mixings: 

BR{N^^m) = -^^-^^^^ , £ = e,//,r . (24) 

We now consider in turn each decay into a specific lepton flavor £. In fig. [1] we 
accordingly display the partial lifetime of A^i , = 1 /F^^ , as a function of its mass and 
for selected values of Ci, its mixing with the light lepton £. We infer that a necessary 
(not sufficient) condition in oder for the heavy neutrino to be stable on cosmological time 
scales is that all the C^'s must be smaller than 10~^^, so that is bigger than the age 
of the universe. 

The primary £^ and W^, being the decay products of a particle nearly at rest, are 
practically monochromatic. In the approximation rrii <^ My/^ we have 

where p and E denote particle's momentum and energy, respectively. 

The decaying heavy neutrinos A^i are unpolarized, so the lepton £ and the W bosons 
are emitted isotropically. Focusing on a single decay, it turns out that f"*" (£~) is emitted 
preferentially in the same (opposite) direction as the A^^i spin. If the decaying A^^i has a 
TeV-scale mass, f"*" is approximately an helicity eigenstate with eigenvalue +1 (— 1). 
The associated W boson can then be considered to be approximately an eigenstate of 
the component of the spin along its direction of motion with eigenvalue, that is to be 
longitudinally polarized. 
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FIG. 2: Example of charge conjugated Feynman diagrams for a cascade decay of A^i. Clearly, this 
is just a suggestive pictorial representation because me > mi' > mg" is realized only for r, fj,, e 
respectively. 



B. Electron and positron energy spectra 

Clearly, if £^ is an unstable lepton like /u^ or r^, it decays with a lifetime much shorter 
than tni, producing and neutrinos. Also the associated boson quickly decays via 
electroweak interactions. Since we are interested in a TeV-scale Ni, its decay products are 
highly relativistic particles and we can safely work in the narrow- width approximation, 
i.e. we treat intermediate particles in fig. [2] as if they were on-shell. Accordingly, the 
number of per unit energy obtained from the decay of A'^i can be written as a sum of 
the contributions from the primary lepton and W boson decay chains: 



where the factor 1 /2 has been introduced in order to fix the normalization of the separate 
contributions to one. 

1. The lepton £ decay chain 

Let analyse first the electrons and positrons coming from the primary £ decay chain. 
£ = e In this monochromatic primary decay products of A^i, see eq. (j25p : 



dWe ~ 2 i dWe ^ dWe 



) 



(26) 




6{E, 



(1 




)) 



(27) 



dEe 



2 
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This is represented in the top right panel of fig. [3] as a vertical line. 



/I In this case, e are produced as secondary decay products from the primary rel- 
ativistic fi^ via its decay into e^Ve^n- They energy spectrum is thus continuous 
between < Eg < E^, as expected for the 3-body decay of a relativistic particle. 
Let take the z-axis to be parallel to the muon momentum. As already stressed, the 
spin of /x+ (^~) is parallel (antiparallel) to z. This means that the /i+ is po- 
larized in its rest frame, having spin parallel (antiparallel) to z, so that the angular 
distribution of the outgoing e"^ (e~) is enhanced in the same (opposite) direction as 
the muon spin, namely the z direction: 



d cos Of. 



_|i(i + _cos^;±) 



(28) 



where ^g± is the angle between z and the momentum of and quantities evaluated 
in the muon rest frame are indicated with a tilde. A relativistic produced in the 
rest frame of the /x^ with energy £'g± , is found in the rest frame of A^i to have an 
energy £'g± = £'e±7(l + cos ^g±), with 7 = E^/m^. As a result, the emitted with 
9f,± < tt/2 are slightly harder than for an unpolarized fi^ and, assuming E'^ ^ m^, 
their number per unit energy is given by 



dN^^'^\Ef, 
dEf.± 



BR{^i^ ^e^vev^) f±{E^±,E,±,±l) , Q < E,± < E^± , (29) 



where E^± is the same as in eq. (j25]) . the BR is close to 100% and the function f± 
is defined as 



f^{E,,Ef,P,) = ^ 



+ « ( §- 



(30) 

where Ei and Ej are the energy of the decaying and final lepton respectively and Pi 
is the polarization of the decaying lepton (see e.g. [39(]). In our case P^± = ±1 and 



since f^{Ei,Ef,l) = f_{Ei,Ef,—l) = f{Ei,Ef), it can be directly checked that 
the energy distributions are equal for electrons and positrons, as expected. We then 
omit the apex it and write explicitly 



dAf(^)(^g 
dE~ 



f{E^,Ee) , f{Ei,Ef) 



4 



1 



(El 

\E, 



(31) 



The energy distribution dN^^\E(.) / dEf, is represented with a dotted curve in the 
bottom left panel of figEl for M^r^ = 1, 2, 3, 5 TeV from top to bottom. 



_r In this case, there are both leptonic and hadronic decay modes, with BR of about 
37% and 63% respectively. 
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As for the leptonic decays, secondary are produced via — t- e^Ve^r with an 
energy distribution similar to eq. pip but suppressed because of the different BR 

TF^^ = ^^(^ ^ ez^ez.-) fiEr,Ee) , 0<Ee<Er (32) 

where E-,- is the same as in eq. (j25p and the BR is about 17.8%. The primary 
decays also into ii^f^u-j- with a similar BR. The latter ji^ has helicity close to 
±1 and fully decays into e^UeU^, giving an additional contribution to the energy 
distribution 



BR{t ^ fliy,,Ur) / dE^f{Er,E^)f{E^,,Ee) (33) 



where Er is as in eq. (l25p and < Ee < Er- The sum of the energy distributions 
resulting from leptonic decays is represented with a dash-dotted curve in the 
bottom right panel of figlS for Mtvj = 1,2,3,5 TeV going from top to bottom. 
Notice that, because of the additional contribution from eq. (|33p . the dash-dotted 
curve is softer than the dotted one in the plot at left. 

Also the hadronic decay modes can produce , which are generically softer than 
those produced in its leptonic decay modes. It is not difficult to give an estimate of 
their energy distribution. The decay modes into one or two pions are ir^Vr and 
P^Vt TT^TT^Ur, with BR~ 11% and 25.5% respectively. The vr^ then fully decays 
into fi^i'fj,. The decay modes into three or more pions have BR of about 26% 
and can be neglected for the sake of the present analysis because the produced 
are very soft. 

Let study first the decay — t- tt^P,- (r~ — ?• 'K~Vr). Since the r"*" (r") has helicity 
close to +1 (—1), in its rest frame it is approximately polarized with spin parallel 
(antiparallel) to z, so that the pion momentum is mostly antiparallel to z 

dt V 

— (1-COS0,) , (34) 



d cos 9-,^ 2 

where 9t^ is the angle between z and the pion momentum in the pion rest frame. 
The pion mostly goes backwards with respect the the tau momentum and therefore 
it is softer than for an unpolarized tau decay. Indeed in the A^i rest frame 



dET^ Pr \ Er 

The vr"*" (vr") decays isotropically in p'^i'^ (^~z?^), so that in the A'^i rest frame the 
muon has simply a flat energy distribution equal to l/pn ~ IZ-E'tt- The (/^~) has 
approximately helicity -|-1 (-1) and fully decays into e^Ve'i^n (e~S'e'^/x)- The energy 
distribution of (e~) is finally: 

^jY(r/7r/M/e)(^^) y-S. fE. I 



dE, 



[ ^ dE /^^}^-> [ \E,-f{E,,E,) . (36) 

JEe dE^ Je, Ptt 
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Secondly, we consider the decay r"*" p'^i^r {t~ — p~Vt)- The p having spin 1, it 
can take several polarization states and its angular distribution is less sensitive to 
the tau polarization 

dt t 

" '1- 0.46 cos 0p) , (37) 



d cos 9p 2 

but still the p preferentially goes backwards with respect to the tau momentum and 
in the A'^i rest frame 

^^ = i?i?(r^p..)A(0.73 -0.46^) . (38) 

The p^ is mildly polarized and it decays nearly isotropically in vr^vr'^, with vr^ — )• 
ji^Vfi- For the sake of the present analysis the electron energy distribution can be 
approximated by 

^ = / dE,"-^^ dE^^ dE,^fiE,,E,) . (39) 
dEe Je, dEp Je, Pp Je, Ptt 

As already remarked, the remaining decays of the involve three or more pions in 
the final state: the eventually produced are even softer and can be neglected in 
the present analysis. 

The sum of the energy distributions resulting from leptonic and hadronic 
decays is represented with a dotted curve in the bottom right panel of figlS for 
MjVi = 1, 2, 3, 5 TeV going from top to bottom. Clearly, the inclusion of the hadronic 
tau decay modes has the effect of softening the spectrum. 



2. The W boson decay chain 

Also the primary and approximately linearly polarized W boson further decays. The 
number of per unit energy produced in its decay chain has to be added to the ones 
studied above. The W leptonic decay modes into pairs ii'i have a BR of about 33%; the 
remaining 67% is represented by hadronic decay modes. Let now identify direction of the 
W momentum with the z axis. 

As for the leptonic decay, in the W rest frame (were the polarization vector is ap- 
proximately e = iz) the secondary i (see figH]) is emitted with angular distribution given 
by: 

- — ;— (1 - cos 9e) , (40) 



dcos Of 



where Oi is the angle between the z axis and the i momentum. The lepton is preferentially 
emitted perpendicularly to z axis in the W rest frame. In the A'^i rest frame it then displays 
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an energy spectrum peaked at E\y/2: 

^^^^^'^ -BRiW^ei.e)—^(l-^] , 0<Ee<Ew (41) 



dEi pw Ew V 

where pw and Ew are given in eq. (j25p and we recall that the BR is roughly 11% for all 
lepton flavors. At this point we have to sum over the flavors i = e, fi, r. When i = e, its 
energy spectrum can be directly read out from eq. (|4ip above: 

^Ni-">(E,) _ ^^(^ ^ A # f 1 - , , (42) 



dEf, py/ E\Y \ Ey/ 

When i = n, the contribution to the e energy spectrum is given by 

''"'"''"'"<^'' - . (43) 



dEf. ./k, dE^ 



When i = T, there is a similar but suppressed contribution 

- BR[t ^ eVeVr) / dEr — f{Er,Ee) . (44) 



dEe JEe 

Since the r also decays also into ^v^Vr, there is another small contribution from 

,^ ^ = BR{t ^ ^iv.vr) / dE,^J^^p± / dE,f{E,,E,)f{E„E,) . 

(45) 

Summarizing, for each primary W^, the total number of produced in its leptonic decay 
modes is approximately (11 + 11 + 2 x 0.17 x 11)% 1/4. 

The hadronic decays of the W amounts to about 67%. They also can produce leptons, 
with a soft spectrum that is negligible for the sake of the present analysis. 

In the top left panel of figl3] we show the number of per unit energy produced by 
summing all the leptonic decays of the (the dotted line is the sole £ = e contribution) , 
which has to be summed to the one coming from the primary and its eventual decays, 
as we are now going to discuss. 



3. Results 

For each lepton flavor, we now sum the energy distributions of from the £ and W 
boson chains. 

For £^ = e^, this sum is easily done since the primary are monochromatic, see the 
top right panel of figlSl 

For £^ = /i^, the energy spectrum of the secondary coming from its decay is soft. 
This is shown by the dotted curve in the left bottom panel of figlSl The solid lines are 
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3' lO* 

Eg Eg 
FIG. 3: Number of per unit energy for selected values of Mni in GcV. Top left: the solid 
lines are the total contribution from the W boson leptonic decay chain; the dotted lines are the 
contribution from W — > ev^. only. Top right: solid lines are the sum of the monochromatic e and 
the W chain. Bottom left: the solid lines are the sum of the ^ and W chains; the dotted ones 
are the contribution of the /i chain only. Bottom right: solid lines are the sum of the r and W 
chains; the dash-dotted lines are the contributions from the leptonic r decays; the dotted lines the 
contribution from both leptonic and hadronic r decays. 



obtained by adding the contribution of the coming from the W chain. Since the W 
chain provides approximately 25% of the total number of produced in this decay of A'^i , 
the contribution is small, especially for the hardest e^. 

For £^ = T^, the energy spectrum of the coming from its leptonic decay is also 
soft, but suppressed because of the BR of about 34%, as can be seen from the dash-dotted 
curves in the right bottom panel of figEl The dotted line are obtained by adding the 
contribution from the r hadronic decay. Finally, the solid lines are obtained by adding 
also the contribution from the W chain, which is relevant at the highest energies and 
hardens the spectrum close to its end point. 
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C. Contributions to CR's 



Having obtained explicit expressions for the number of electrons and positrons per unit 
energy produced via each channel A''i — )• £W and its subsequent decays, it is straightfor- 
ward to substitute them into eg . (1201) and to calculate the flux observed at Earth 

<l>{E,)=Y,cP\E,) , = , (46) 

e 

where c is the speed of light. 

Clearly, because of the propagation in the galaxy, the energy spectrum of the electrons 
and positrons gets modified according to the diffusion equation. The solid lines in the 
three panels of figjH display the combination E^cl)^{Ee) for each lepton flavour, by assuming 
CiX^Jj"^ = 10^^"^ and for four selected values of Mjv^ in GeV. We recall that the fluxes are 
proportional to C|. It is thus not difficult to obtain the total flux (j){Ei.) corresponding 
to any possible set of Ce,Ci^i and Cr- The dotted line are obtained by neglecting the 
contribution of electrons and positrons originated from the W boson chain. The spectrum 
turns out to be hard for i = e, while it is softer for i = and even more for i = t. 

Notice that energy losses due to the diffusion in the galaxy are more relevant at low 
rather than at high E^,. This can be realized by comparing, for instance, the plots for the 
muon case: the value of the flux increases with Mjv^ in the whole energy range, while this 
was not the case for the number of electrons and positrons per unit energy at the source, 
see figlSj The same applies for the tau case, for which the energy spectrum at the source is 
even softer. This justifies the fact that in our calculation we neglected the softer electrons 
and positrons coming from the hadronic r decays with more than three pions, as well as 
those coming from the hadronic decays of the W. Fig|3]show indeed that the shape of the 
fluxes after propagation is essentially controlled by the behaviour of dN/dEe close to its 
end point {i.e. for the highest possible values of Ef.). For the muon case, this behaviour 
is a direct reflection of eq. ([3T]) : for the tau, of the sum of eqs. ([32]) . ([33|) and P2]) . 



The same argument can be applied to electroweak radiation effects |25|, l28|. Since the 
mass scale of the Majorana neutrino A'^i is much bigger than the electroweak scale, soft 
electroweak vector bosons can be radiated in its decay. Such vector boson give even softer 
electrons and positrons, they can be neglected in the present analysis because the cosmic 
ray excesses observed by PAMELA and Fermi-LAT are comprised in the energy range 
from 20 GeV to 1 TeV. 



IV. COMPARISON BETWEEN THE MODEL AND THE DATA BY PAMELA 
AND FERMI-LAT 

We now come to the comparison between our fourth lepton family model with a decay- 
ing TeV-scale Majorana neutrino and the experimental data on CRs recently collected 
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e 





by PAMELA and Fermi-LAT. We will first describe an extremely fast method, based on 
the Sum Rules (SR) introduced in ref. [261] . This method allows for a fast scrutiny of the 
viability of any model (we considered in particular the case of models predicting an equal 
amount of fluxes for electrons and positrons, but the method could of course be general- 
ized). Then, we go through the longer path of the direct comparison, also to show the 
reliability of the SR method. 



A. Sum Rules method 



A fast method to check the viability of our model is to make use of the results of 



ref. 26|, that we briefly summarize here. For definiteness and also for an easy compari- 
son with the literature, let consider the Moskalenko Strong (MS) 35|, l3g| fluxes B^{Ee), 
where it refer to positrons and electrons respectively, to model the contribution of astro- 
physical sources (one could of course apply the same reasoning for any other astrophys- 
ical background model), leaving the normalization of the fluxes A'^^ as a free parameter: 

ckl{Ee)=NBB^{Ee). 
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We assume that an unknown source of CR is also present, which gives equal fluxes 
for electrons and positrons, (j)^{Ee) = (fP_{Ef>) = ((P{E(.). Hence the total fluxes are 

MEe) = (j)"^ {Ee) + (t>l{Ee) , 0_ (S^) = .^^ (^e) + (£^e) • (47) 

In ref ■ j26l| we showed how to combine the PAMELA and Fermi-LAT data in order to 
extrapolate the flux of from this unknown source. We re-display this flux (p^ (Eg) in 
fig- El by means of the shaded region. The inner and outer bands have been obtained 
by combining respectively the la and 2a error bands of PAMELA and Fermi-LAT, while 
considering all possible associated values for Nb by means of the SR. 

It is now straightforward to compare the experimentally extrapolated flux cjp {Ee) with 
our particular model, belonging indeed to the category of models which predict equal 
fluxes for the electrons and positrons. All we have to do is to superimpose the (solid) 
curves of the fluxes already shown in figlH This is done in fig. [5l where the C^'s have 
been fixed so to cross the shaded region at the lowest energies. We recall that the fiuxes 
scale as C|, so that it is not difficult to obtain the value of the fiuxes for any other value 
of C^. The dashed curves have been obtained by neglecting polarization effects in both I 
and W decay chains. 

We can see that the electron and positron fluxes obtained via the mixing of A^i and the 
e-flavor charged lepton have not the right spectral index to fit the combined experimental 
data. In the case of mixing with the muon, the electron and positron fiuxes follow the 
shape of the data quite well for values of A^^i close to 2 — 3 TeV. For the mixing with tau, 
the agreement is even more impressive. For Mtv^ = 3 TeV, the curves have been obtained 
by setting Cgx'J.l = IQ-^^-oa, c^x^J,^ = IQ-^^-H and Crx'J.l = IQ-^^-OS. 

This means that, if we want to interpret the PAMELA and Fermi-LAT excesses as 
being due to a decaying heavy Majorana neutrino, we need M^^ 3 TeV, CrX^^"^ ^ 10"^'^ 
and/or C^j^x^^"^ w 10"^"^ while CeX^^"^ < 10~^^. Hence the hierarchy Of, < O^^r among the 
mixing angles with the light leptons is needed. 

The universal coupling case_^e = G^i = Gr is thus excluded, as already pointed out by 



the numerical analysis of ref.[2j], where a generic heavy fermion decaying into £W was 
considered (the nature of the vertex interaction was not specified). We postpone to the 
next section a more complete comparison with previous analysis. 



B. Direct method 



We now come to the direct comparison between our fourth lepton family model with a 
decaying TeV-scale Majorana neutrino, considering separately the PAMELA and Fermi- 
LAT experimental data on CRs. 

The fiuxes of arising from A^i decays that we come to calculate have to be summed 
to those arising from astrophysical sources. As before, we consider the popular MS model 
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£, [GeV| E, [GeV] 

FIG. 5: Comparison between the prediction of our model for selected values of Mni in GeV and 
the combined data from PAMELA and Fermi-LAT 26[. The inner and outer bands have been 



obtained by combining respectively the Icr and 2a error bands of PAMELA and Fermi-LAT (2010 
releases), together with the MS model for the astrophysical background. For M^-^ = 3 TeV, the 



curves have been obtained by setting CeX 



1/2 

Ni 



10 ^^-^^j C^Xpf^ 



1/2 



10-2711 and Crx]l^ 



10 



-27.05 



Solid (dashed) lines display the results of including (neglecting) polarization effects. 



35l . l36l | for such astrophysical background and recall the explicit expressions for B^{E): 



1 + 650£;2-3 + 1500£;4.2 ' 

1 + 11^0-9 + 3.2^2.15 ' 

1 + UOE^-^ + 600-E;2-9 + 580E^-^ ' 



(48) 
(49) 
(50) 

(51) 



where E is measured in GeV and the Bs in GeV~^cm-2sec-^sr-^ units. To be specific, 
in the present analysis we fix the value A'^^ = 0.638, obtained by means of the sum rule 
method proposed in ref. 26l |. 

We start by considering the positron fraction and display it in the left panels of fig. [6] 
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for i = e,iJ,,T respectively and for selected values of the A^i mass, Mjy-^ = 1,2,3,5 TeV. 
Solid (dashed) curves have been obtained by including (neglecting) polarization effects. 
The panels also show the PAMELA la data points (2010 release) and the MS model 
prediction. The mixing angle Ci has been chosen so that our model curves cross the 
central value of the PAMELA positron fraction data point associated to the 35 GeV 
energy bin, 4>e+/{4>e~ + 4'e+) = 0.733. Notice that plays the role of shifting up and 
down the curves with respect to MS. The shape of the curves depends on the value of 
Mat^, since they drop at about MjVi/2. The dropping is sharp for £ = e, soft for and 
even softer for r. This pattern clearly reflects the shape of the fluxes displayed in figlH 
It can be realized that the PAMELA datapoints are nicely fitted for any value of M/v^, 
provided it is bigger than about 200 GeV (8|. 

This does not happen for the Fermi-LAT data on the total electron and positron flux, 
as shown in the right panels of fig. [6] for £ = e,fi,T respectively. The couplings Ci have 
been chosen to be the same as the corresponding ones in the PAMELA plots at the left. 
Again, solid (dashed) curves have been obtained by including (neglecting) polarization 
effects. It is evident that the Fermi-LAT data points cannot be fitted for i = e, since the 
curves raise too much. For (. = t the curves have a shape that is nicely compatible with 
the data for values of M^Vj close to 2 — 3 TeV, while for i = fi there is some tension because 
of an excess around 500 — 1000 GeV. For Mtvj = 3 TeV, the values of the combination 

1/2 

logiQ{C£X^_^ ) used in the plot are —27.105, —27.112, —27.031 respectively for e,iJ,,T. This 
means that, if we want to interpret the PAMELA and Fermi-LAT excesses as being due 
to a decaying heavy Majorana neutrino, we need 6e < dfi,T, as concluded in the previous 
section. 

As already stressed, the curve for r is slightly softer than for fi because the its energy 
spectrum is slightly softer, see figl3] and HI The excess of the curves at 2 — 3 TeV 
could be easily corrected for by slightly reducing the coupling C^. This is demonstrated 
in figH where we show the best fit values of Ce for i = ^ and r and for M^r^ = 3 TeV. 
We can see that the r fits astonishingly well both PAMELA and Fermi-LAT data with 
x^Jj'^Cr = 10"^^''''^ — lO"^''''^^. The fit for the muon with a slightly suppressed coupling, 
x^jlj'^C^ = 10"^'^-^^ — lO"^'''^'^, is now good. Clearly, we were able to conclude all this also 
by means of the SR method discussed previously. However, the direct analysis allows for 
a better comparison with the results of ref. 2J], as we now turn to discuss. Polarization 
effects have been included in the upper panels, while they have been neglected in the lower 
panels. The inclusion of polarization slightly hardens the total flux. 

Ref. 2M concluded that the tW channel does not fit the Fermi-LAT data because of a 



too much flat curve; their best flt was obtained in the case of the channel. Indeed, 
in their plots the fluxes appear softer for both and tW as compared to figlZl even 
neglecting polarization effects; the latter is in particular so soft that the agreement with the 
Fermi-LAT experimental data is lost. In ref. [Ul the tW mode was considered to provide 
a good fit of the data, but worse than the jiW mode. We suspect that the discrepancy 
with previous numerical analysis could either be due to a vertex interaction different from 
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FIG. 6: Comparison with the 2010 releases of PAMELA (left) and Fermi-LAT (right), considering 
separately the mixing between iVi and each SM lepton t = e,fi,T. The value of the Ci fits the 
central value of the PAMELA 35 GeV data point. The solid lines correspond to — 1, 2, 3, 5 
TeV from bottom to top respectively. The shaded region emphasizes the energy range where 
solar modulation effects cannot be neglected. The bottom black line displays the MS estimate for 

1/2 

the astrophysical backgrounds. For Mni — 3 TeV, the values used in the plot for \ogiQ{CeXj^^ ) 
are —27.105, —27.112, —27.031 respectively for e,iJ,,T. Solid (dashed) lines display the results of 
including (neglecting) polarization effects. 
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WITHOUT POLARIZATION 



WITHOUT POLARIZATION 





E, [GeV| E, [GeV| 

FIG. 7: Our best fits of PAMELA (left) and Fermi-LAT (right) data. Tlie mass of Ni lias 
been chosen to be 3 TeV. In the upper (lower) panels polarization effects have been included 
(neglected). For A^i tW wc used x]lfCr = IO^^t.ot _ 10-27.02^ ^YyWc for iVi fiW we used 

X 



'-jvf = 10 ^^'^^ ^ 27.20 jj-^ ^Yic upper panels polarizations effects are included, in the lower 
panels they are neglected. 



eq-dH) (for instance including all chiralities) or to their numerical code. Their analysis is 
based on PYTHIA, that treats leptons and vector bosons as unpolarized. However, as 
diplayed in figlTl the inclusion of polarization effects cannot alone explain the discrepancy. 

The hardness of the upper right panel spectra for fi and r follows directly from the 
fact that the CC vertex involves L-handed (R-handed) leptons (antileptons). If we had 
considered a non chiral vertex or if we had neglected polarization effects in our calculation, 
the curves in figJT] would have been softer, as demonstrated in the lower right panel. This 
can be easily understood for the /i case, since the contribution from the W decay chain 
is negligible and the shape of the flux in figiT] is essentially controlled by the energy 
distribution of the produced in the /j,^ decay, eq. (j3T]) . For instance, for a relativistic 
L-handed /U~, the highest energy e~ produced in its decay are slightly harder than in the 
unpolarized case because they are preferentially produced with momentum parallel to the 
muon one, eq. (j28p . At the contrary, for relativistic R-handed fi~ , the highest energy 
electrons preferentially escape with opposite momentum and are thus softer. 
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Ej, [GeV] E„ [GeVI 

FIG. 8: Antiproton flux (left) and antiproton/proton ratio (right) for the MAX, MED, MIN 
propagation models, = 3 TeV and x^j^Ci = 10^^^'^. The dashed curves in the left panel 

show the primary antiprotons for MAX, MED, MIN models. The lower shaded region represents 
the secondary antiprotons only, according to ref . [40I1 . The upper three curves correspond to the 
sum of the primary and secondary antiprotons. The (lower) PAMELA data have smaller error 
bars with respect to the (upper) CAPRICE data Q. 

V. COSMIC RAY ANTIPROTONS 

In this section we give an estimate for the cosmic ray antiproton flux resulting from 
the heavy Majorana neutrino decay and compare it with current data. 

Protons and antiprotons are generated via the hadronic decay of the primary W boson, 
with BR of about 67%. Their energy spectrum is determined by the fragmentation and 
hadronization processes. In this case an analytic approach is not suitable and we rather 



adopt the numerical recipies provided in ref. j25l |. In particular, the antiproton (proton) 
flux obtained from a 1.5 TeV dark matter annihilating into W^W~ is twice the antiproton 
(proton) flux in our model. 

The propagation for antiprotons through the galaxy is described by a diffusion equation 



whose solution can be cast in a factorized form as discussed e.g. in [25[. In this case the 
astrophysical uncertainties associated to the dark matter profile and to the propagation 
parameters are large, about one order of magnitude. We display the antiproton flux in the 
left panel of fig. Elfor the MAX, MED, MIN propagation models and for Mn^ =3 TeV and 
x]I^Ce = 10-27-2. The dashed curves show these primary antiprotons for MAX, MED, 
MIN models; the lower shaded region represents the flux of the secondary antiprotons 



according to ref. 40(]; the upper three curves correspond to the sum of the primary and 



secondary antiprotons. The (lower) PAMELA data [6|] have smaller error bars with respect 



to the (upper) CAPRICE data [27[. It turns out that only the MIN propagation parameter 
set is compatible with the data, the MED one is barely compatible, while the MAX seems 
disfavored. 

The antiproton/proton ratio is studied in the right panel of figlHl For the proton 
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flux, we consistently adopt the parameterization of ref. 40|] with spectral index equal to 
—2.72. This parameterization is valid for energies higher than about 10 GeV. Again, 
the (lower) PAMELA data points 0] are more precise than the (upper) CAPRICE ones 



27|. The plot shows that even the MIN model displays tension with the PAMELA data. 



Future measurements confirming the PAMELA results could be able to rule out the heavy 
decaying neutrino model studied in this work. 



VI. CONCLUSIONS 



In this work we have asked whether the electron and positron CR excesses observed by 
PAMELA and Fermi-LAT can be produced by a decaying fourth lepton family Majorana 
neutrino A'^i. Several more or less natural models of Majorana neutrino dark matter have 
been put forward in the literature in addition to the ones cited in this work. A review of 
different models is We note, however, that our model of Majorana neutrino emerges 
naturally at the electroweak scale given that the new generation of leptons is needed 
to resolve the Witten anomal y a ssociated to the specific technicolor extension known as 
Minimal Walking Technicolor 16l4l8l| . We have already investigated the collider properties 
of our extension in [19|]. There we have also discussed and compared with other models of 



fourth family of leptons at the electroweak scale. 

This heavy neutrino mixes only with the SM charged leptons via CC interactions, so 
its primary decay products are £W, i = e, fi, r, with BR proportional to the combination 
C|/(C'e + + C*^)' where Ci represents the strength of the mixing between N and the 
SM charged lepton i. If the mixing angles are sufficiently small, A'^i is nearly stable on 
cosmological time scales and constitutes a fraction x of the today dark matter density. 

We analytically calculated the energy spectrum of the electrons and positrons produced 
via the decay channel iW and its subsequent decays, including polarization effects. The 
shape of the energy spectrum is different for i = e, fj,, r, as can be seen in figlSl 

After propagationfhere we adopted the approximation of 2^ for the Green function 
of the MED-model [Sg]) the contribution to the total flux from each decay channels can 



be straightforwardly compared with the PAMELA and Fermi-LAT experimental data 
combined according to the Sum Rules method [261 ]. see figEl We also carried out the more 
standard method of the separate comparison between the PAMELA and Fermi-LAT data 
in order to show the reliability of the SR method. 

Our results are that the model provide a plausible explanation of the observed CRs 
excesses, when the BRs of A^i — t- tW and A'^i — t- fj,W are the dominant ones, as it is clear 
from fig. [71 No acceptable agreement with the data can be achieved for A'^i — t- eW. This 
means that A'^i should have a hierarchical pattern of mixing angles with the SM charged 
leptons, with Cg < C^,r- In particular, the tW decay mode alone seems to best fit the 
data. In this case, the best fit occurs for M^^ 3 TeV and Afi Cr ~ 10"^"^. For the tW 
decay our results differ from the ones obtained in ref. [ill, by using numerical codes. 
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We found that our inclusion of polarization effects cannot alone explain this discrepancy. 

We have demonstrated that the positron and electron CR excesses possibly present 
in the PAMELA [45| and Fermi-LAT data, can be simply accounted for in our fourth 
generation heavy Majorana neutrino model. Clearly, the contraints from gamma pho- 
tons [4^, |43[ would deserve a dedicated analysis, but we do not expect to find tension 



between our model and these data. We instead estimated the primary contribution to 
the antiproton flux and the antiproton to proton ratio in our model. Despite the uncer- 
tainty in the propagation parameters, it seems that even the MIN model shows tension 
with the PAMELA data fl. The AMS-02 space station experiment will hopefully provide 
additional relevant informations 44 1. 
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Appendix A: Full mixing 

The mixing between the neutrinos of flavor (" and the three SM neutrinos is described 
by the 5x5 mass matrix A4: 
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Barring unnatural tunings and up to corrections to its mixings of 0{eV/MN^^N2), the 
unitary matrix V is 
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where 
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diagonalises the lower 2x2 sector of the 5x5 mass matrix Ai , namely 



/ eff 



VMV^ 



m 



3x3 











\ 



V 



Mn^ 
J 



(A4) 



where the elements of Tn'^{/;^ are naturally 0{eV) and 
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At this stage, the matrix in eq. (|A4p can be fully diagonalised with a further unitary matrix, 
which can be identified with the MNS mixing matrix, acting only on the upper 3x3 block: 



U^ml{{u* = diagi 



mi, 1712,1713) 



(A6) 



Notice that in the limit of small mixing angles between the heavy and light neutrinos, the 
mixing matrix V becomes (at first order in se) 



V 



icsp ss, 
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(A7) 



We have thus demonstrated that each light flavor couples to the heavy neutrinos as in 
eq.([l2]). 



28 



[1] O. Adriani et al. [PAMELA Collaboration], "An anomalous positron abundance in cosmic rays 
with energies 1.5-100 GeV, Nature 458 (2009) 607 [arXiv:0810.4995 i [astro-ph]]. O. Adriani et 
al. [PAMELA Collaboration] , A statistical procedure for the identification of positrons in the 
PAMELA experiment, Astropart. Phys. 34 (2010) 1 |arXiv:100l!3522l [astro-ph.HE]]. 

[2] A. A. Abdo et al. [Fermi LAT Collaboration], Measurement of the Cosmic Ray e+ plus e- 
spectrum from 20 GeV to 1 TeV with the Fermi Large Area Telescope, Phys. Rev. Lett. 102 
(2009) 181101 jarXiv:0905.0025l [astro-ph.HE]]. M. Ackermann et al. [Fermi LAT Collabora- 
tion], Fermi LAT observations of cosmic-ray electrons from 7 GeV to 1 TeV, Phys. Rev. D 
82 (2010) 092004 |arXiv:1008.3999l [astro-ph.HE]]. 

[3] J. Chang et al. [ATIC Collaboration], An Excess Of Gosmic Ray Electrons At Energies Of 
300-800 Gev, Natm-e 456 (2008) 362. 

[4] S. Torii et al. [PPB-BETS Collaboration], High-energy electron observations by PPB-BETS 
flight in Antarctica. larXiv:0809. 07601 [astro-ph]. 

[5] F. Aharonian et al. [H.E.S.S. Collaboration], The energy spectrum of cosmic-ray electrons at 
TeV energies, Phys. Rev. Lett. 101 (2008) 261104 ja:rXiv:0811.3894l [astro-ph]] . F. Aharonian 
et al. [H.E.S.S. Collaboration], Probing the ATIG peak in the cosmic-ray electron spectrum 
with H.E.S.S, Astron. Astrophys. 508 (2009) 561 iarXiv:0905.0105l [astro-ph.HE]] . 

[6] O. Adriani et al. [PAMELA Collaboration], PAMELA results on the cosmic-ray antipro- 
ton flux from 60 MeV to 180 GeV in kinetic energy, Phys. Rev. Lett. 105 (2010) 121101 
[arXiv: 1007.08211 [astro-ph.HE]]. See also: O. Adriani et al., A new measurement of the 
antiproton-to-proton flux ratio up to 100 GeV in the cosmic radiation, Phys. Rev. Lett. 102 
(2009) 051101 ; arXiv:0810:4994l [astro-ph]] . 

[7] Y. -Z. Fan, B. Zhang, J. Chang, Exec sses in the Gosmic Ray Spectrum and Possible 
Interpretations, Int. J. Mod. Phys. D19 (2010) 2011-2058. [arXiv: 1008.46461 [astro-ph.HE]] . 

[8] E. Nardi, F. Sannino and A. Strumia, Decaying Dark Matter can explain the electron /positron 
excesses, JCAP 0901 (2009) 043 jarXiv:0811.4153l [hep-ph]]. 

[9] A. Ibarra, D. Tran and C. Weniger, Detecting Gamma-Ray Anisotropics from Decaying Dark 
Matter: Prospects for Fermi LAT, Phys. Rev. D 81 (2010) 023529 iarXiv:0909.3514l [hep-ph]] . 
[10] M. Cirelli, P. Panel and P. D. Serpico, Diffuse gamma ray constraints on annihilating or 
decaying Dark Matter after Fermi, Nucl. Phys. B 840 (2010) 284 |arXiv:0912!0663 [astro- 
ph.CO]]. 

[11] P. Meade, M. Papucci, A. Strumia and T. Volansky, Dark Matter Interpretations of the 
Electron/Positron Excesses after FERMI, Nucl. Phys. B 831 (2010) 178 [ arXiv:0905.0480l 



29 



[hep-ph]]. 

[12] M. Papucci and A. Strumia, Robust implications on Dark Matter from the first FERMI sky 



gamma map, JCAP 1003 (2010) 014 |arXiv:0912.0742 [hcp-ph]]. 
[13] G. Hutsi, A. Hcktor and M. Raidal, Implications of the Fermi-LAT diffuse gamma-ray mea- 
surements on annihilating or decaying Dark Matter, JCAP 1007 (2010) 008 [ arXiv: 1004^20361 
[astro-ph.HE]]. 

[14] C. Boehni, T. Delahaye and J. Silk, Can the morphology of gamma-ray emission distinguish 
annihilating from decaying dark matter?, Phys. Rev. Lett. 105 (2010) 221301 [ arXiv:1003. 12251 
[astro-ph.GA]]. 

[15] S. Palomares-Ruiz and J. M. Siegal-Gaskins, Annihilation vs. Decay: Constraining Dark 
Matter Properties from a Gamma-Ray Detection in Dwarf Calaxies, larXiv: 1012.2335 ^ [astro- 
ph.CO]. 

[16] F. Sannino and K. Tuominen, Technionentifold, Phys. Rev. D 71, 051901 (2005) 
|arXiv:hep-ph/040520^ . 

[17] D. D. Dietrich, F. Sannino and K. Tuominen, "Light composite Higgs from higher represen- 
tations versus electroweak precision measurements: Predictions for CERN LHC," Phys. Rev. 

D 72, 055001 (2005) |arXiv:hep-ph/0505059| . 
[18] R. Foadi, M. T. Frandsen, T. A. Ryttov and F. Sannino, "Minimal Walking Technicolor: Set 

Up for Colhdcr Physics," Phys. Rev. D 76, 055005 (2007) arXiv:0706. 16961 [hep-ph]]. 
[19] M. T. Frandsen, I. Masina and F. Sannino, Fourth Lepton Family is Natural in Technicolor, 

Phys. Rev. D 81, 035010 (2010) jarXiv:0905. 13311 [hep-ph]]. 
[20] O. Antipin, M. Heikinheimo and K. Tuominen, Natural fourth generation of leptons, JHEP 

0910 (2009) 018 ;a rXiv:0905.0622l [hep-ph]]. 
[21] J. R. Andersen et al.. Discovering Technicolor. larXiv:1104. 12551 [hep-ph]. 
[22] A. Ibarra and D. Tran, Antimatter Signatures of Cravitino Dark Matter Decay, JCAP 0807 

(2008) 002 [arXiv:0804.4596 [astro-ph]]. 
[23] A. Ibarra and D. Tran, Decaying Dark Matter and the PAMELA Anomaly, JCAP 0902 (2009) 

021 |arXiv:0811.l"555l [hep-ph]]. 
[24] A. Ibarra, D. Tran and C. Wenigcr, Decaying Dark Matter in Light of the PAMELA and 

Fermi LAT Data, JCAP 1001 (2010) 009 [arXiv:0906.i57ll [hep-ph]]. 
[25] M. Circlh et al., PPPC 4 DM ID: A Poor Particle Physicist Cookbook for Dark Matter Indirect 

Detection, JCAP 1103 (2011) 051 |arXiv:1012.45T5l [hep-ph]]. 
[26] M. T. Frandsen, I. Masina and F. Sannino, Cosmic Sum Rules. larXiv:1011. 00131 [hep-ph]. 
[27] M. Boczio et al. [WiZard/CAPRICE Collaboration], The cosmic-ray anti-proton flux between 

3-CeV and 49-CeV, Astrophys. J. 561 (2001) 787 |arXiv:astro-ph/0103513| . 



30 



[28] P. Ciafaloni, D. ComcUi, A. Riotto, F. Sala, A. Strumia and A. Urbano, Weak Corrections 
are Relevant for Dark Matter Indirect Detection, JCAP 1103 (2011) 019 [ arXiv:1009!0224l 
[hcp-ph]]. 

[29] The Review of Particle Physics, K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 
075021 (2010). 

[30] F. del Aguila, J. A. Aguilar-Saavedra and R. Pittau, Heavy neutrino signals at large hadron 

colliders, JHEP 0710 (2007) 047 |arXiv:hep-ph/0703261] . 
[31] D. Tucker-Smith, N. Weiner, Inelastic dark matter, Phys. Rev. D64 (2001) 043502. 



[hep-ph/0101138]. 



[32] W. -Y. Keung, P. Schwaller, Long Lived Fourth Generation and the Higgs, j arXiv :1103.3765l 
[hep-ph]]. 

[33] J. F. Navarro, C. S. Frenk, S. D. M. White, The Structure of cold dark matter halos, Astrophys. 



J. 462 (1996) 563-575. |astro-ph/9508025] . 

[34] M. Kamionkowski, M. S. Turner, A Distinctive positron feature from heavy WIMP annihila- 
tions in the galactic halo, Phys. Rev. D43 (1991) 1774-1780. 

[35] A. W. Strong, I. V. Moskalenko, Propagation of cosmic-ray nucleons in the galaxy, Astrophys. 
J. 509 (1998) 212-228. [ astro-ph/9807150j . A. W. Strong, I. V. Moskalenko and V. S. Ptuskin, 
Cosmic-ray propagation and interactions in the Calaxy, Ann. Rev. Nucl. Part. Sci. 57 (2007) 
285 |arXiv:astro-ph/0701517] . 

[36] E. A. Baltz, J. Edsjo, Positron propagation and fluxes from neutralino annihilation in the 
halo, Phys. Rev. D59 (1998) 023511. |astro-ph/9808243| . 

[37] J. Hisano, S. Matsumoto, O. Saito, M. Senami, Heavy wino-like neutralino dark matter anni- 
hilation into antiparticles, Phys. Rev. D73 (2006) 055004. |hep-ph/0511118| . 

[38] T. Delahaye, R. Lineros, F. Donato, N. Fornengo, P. Salati, Positrons from dark matter 
annihilation in the galactic halo: Theoretical uncertainties, Phys. Rev. D77 (2008) 063527. 
!arXi v:0712.23T2l [astro-ph]]. 

[39] P. Abreu et al. [DELPHI Collaboration] , Measurements of the tau polarization in ZO decays, 
Z. Phys. C 67 (1995) 183. 

[40] T. Bringmann and P. Salati, The galactic antiproton spectrum at high energies: background 



expectation vs. exotic contributions, Phys. Rev. D 75 (2007) 083006 arXiv:astro-ph /06 125 14| . 
[41] O. Adriani et al. [ PAMELA Collaboration ], "The cosmic-ray electron flux measured by 

the PAMELA experiment between 1 and 625 GeV," Phys. Rev. Lett. 106, 201101 (2011). 

farXi v:1103.2880l [astro-ph.HE]]. 
[42] A. A. Abdo et al.. Observations of Milky Way Dwarf Spheroidal galaxies with the Fermi- 

LAT detector and constraints on Dark Matter models, Astrophys. J. 712 (2010) 147 



31 



[arXiv: 1001.45311 [astro-ph.CO]]. 
[43] A. A. Abdo et al. [The Fermi-LAT collaboration], The Spectrum of the Isotropic Diffuse 
Gamma-Ray Emission Derived From First- Year Fermi Large Area Telescope Data, Phys. Rev. 
Lett. 104 (2010) 101101 |arXiv:1002.3603l [astro-ph.HE]]. 



[44] See: http://ams.cern.ch/ 



[45] More recent data of the PAMELA collaboration [4l| on the negative electron flux has been 
released which, however, given the large uncertainties for energies higher than 100 GeV do 
not affect our results. 



32 



